Figure 1 







E 




o 




w 












o 












X 








LL 


O 



10 20 30 40 
Time after 1993 March 23 [day] 











OJ 




E 




o 




w 




CD 




i_ 


0) 


o 












X 








LL 


O 









t 






c 


f * 


t o 

+ 1 < 


t 


) - 























500 1000 1500 2000 
Time after 1993 March 23 [day] 




AT „ O 

No. a 














- 










38.5/14 









~ r~ 
No. 4. 



0.1/$ 

_L 



2.5 



12 12.5 13 13.5 20.5 




1 




No. 7 | 










12.8/7 
i 





No 


.8 


_ 














17 


2/7 





No 


.9 


_ 








36 


5/8 





44.5 202.5 203 361 361.5 564 564.5 565 



No. 11 












22.0/9 







_ 

No 


12 




_ 


















20 


0/9 







727.5 932.5 933 933.5 1107 1107.5 1108 



No. 14 


















41.9/9 







No. 15 



253.4/9 



_L 





No. 16 . 




♦_ 




0.9/3 . 



1302 1494.5 1495 



1831 1831.5 1832 2024.5 2025 



Figure 3 



(a) o 



o 

LL 

CD 

■H — 1 

o 

-I— • 

CO 




(b) 



o 

"■4—" 

o 

LL 

CD 

■H — 1 

o 

-I— • 

CO 



10 
1 

0.1 
0.01 
0.001 



(c) o 



o 



0) 

-I— • 

o 

-I— ■ 



10 

1 

0.1 
0.01 
0.001 



0.1 1 10 100 1000 

Time-lag [day] 



0.1 1 10 100 1000 

Time-lag [day] 





'"I ! ' ' 


h 4 : 

! _ + .« 


aa '■ 

2 . :*s# : 

x ! ■ 
s * • 5 s:8s ° 






?If I, : 














; ++ 
















: ? si*** 
: s t * j ^ 

;: o ♦ 

...■8 i = 


■ . e 

; : 'e » ! 
* i * * ' i 




: : : + ■ : 
: : o : D : ■ 

: : " t>° \ \ 





0.1 1 10 100 1000 

Time-lag [day] 




0.001 



(e) .1 



o 

LL 

CD 

o 0.01 
3 



10 

1 

0.1 



CO 



0.001 




0.1 1 10 100 1000 

Time-lag [day] 

















1 

i f * 




: o ■ 
















- - -* 

















0.1 1 10 100 1000 

Time-lag [day] 



Figure 4 




Figure 5 



(a) 



o 
c\i 

— CD 

Q ^ 

tt) 



CM 







i i ■ ■ 


i i i i 


• •II 




1 1 


. \ 

\ 

_..N 

\ 


' \ 
\ 














\ 






































. . . . 


. . . . 







100 _ 1000 10000 
1 /Break frequency [day] 



(b) 



^10000 



CD 

-a Q 

CD — 
N CO 

CO O 

o if) 



100 
1 



cd 0-01 
o 

°- 0.0001 




0.0001 0.001 0.01 0.1 1 
Frequency [1/day] 



Figure 6 




log (Frequency [Hz]) 



Mon. Not. R. Astron. Soc. 000, 000-000 (2000) 



Printed 1 February 2008 



(MN MfeK style file vl.4) 



Long-Term Variability of the Low-Luminosity Active 
Galactic Nucleus of M81 



Naoko Iyomoto 1 , Kazuo Makishima 2 

1. The Institute of Space and Astronautical Science, 3-1-1 Yoshinodai, Sagamihara, Kanagawa 229-8510, Japan 

2. Department of Physics, School of Science, The University of Tokyo, 7-3-1 Kongo, Bunkyo-ku, Tokyo 113-0033, Japan 



Accepted 2000 September 22. Received 2000 August 14; in original form 2000 April 28. 



ABSTRACT 

Long-term X-ray variability of the low-luminosity active galactic nucleus of M81 was 
studied, using 16 ASCA observations spanning 5.5 years. The object exhibits a factor 
three variation over the 5.5 years. The source intensity was relatively constant within 
each observation which lasted typically for one day, but intra-day variability by 30% 
was detected on the 15th observation. The power-spectral density (PSD) was estimated 
in "forward" manner, over a frequency range of 10 _8 ' 2 -10 -4 3 Hz (period range of 1/4 
day-5.5 years), by utilizing structure function and extensive Monte-Carlo simulations 
in order to overcome the very sparse and uneven data samplings. When the PSD 
is assumed to be white below a "break frequency" /b and falls off as cx f~ a above 
/b, where / is frequency and a is a positive parameter, the M81 light curve is well 
described with l// b > 800 days and a = 1.4 ± 0.2. 
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1 INTRODUCTION 

A number of low- luminosity AGNs (LLAGNs), with X-ray 
luminosities in the range i() 40_42 erg s _1 , have been discov- 
ered in nearby normal galaxies with ASCA (e.g. Makishima 
et al. 1994; Iyomoto et al. 1996; Ishisaki et al. 1996; Iy- 
omoto et al. 1998; Terashima et al. 1998; Ptak et al. 1999). 
LLAGNs are now considered to be ubiquitous at least in our 
neighbourhood. 

The low luminosities of LLAGNs imply lower black-hole 
masses and/or lower mass accretion rates as compared to 
more luminous AGNs. These alternatives may not be dis- 
tinguished by X-ray spectra, since AGNs exhibit similar X- 
ray spectra over wide range of luminosities, from LLAGNs 
through quasars. On the other hand, X-ray variability is ex- 
pected to provide useful clues to the issue, as a zero-th order 
approximation, because characteristic time scales of the X- 
ray emission is considered to correlate positively with the 
system size and hence the system mass. 

During individual observations with ASCA, which typi- 
cally last for one day, LLAGNs tend to be much less variable 
than the lower-luminosity end of the Seyfert galaxies (e.g. 
Ptak et al. 1998) . We are therefore tempted to speculate that 
LLAGNs have longer characteristic time scales of the X-ray 
variation and hence larger black-hole masses compared to 
Seyferts. In this case, LLAGNs might be "dead quasars" in 
our neighbourhood. To examine this possibility, it is impor- 



tant to determine the power spectral density (PSD) of in- 
tensity variations of LLAGNs on time scales such as dozens 
of days or longer. So far, no such attempts have been made, 
except for a rough estimation by Ishisaki et al. (1996), be- 
cause only few observations with uneven intervals are usually 
available for each LLAGN in X-rays. 

Using ASCA, we investigate long-term variability of an 
LLAGN at the nucleus of M81, which is an Sab galaxy at 
a distance of 3.6 Mpc (Freedman et al. 1994). Its nucleus 
exhibits weak AGN properties in various wave bands. Op- 
tically, it is classified as a LINER (Ho et al. 1997). In ra- 
dio, the spectrum from the nuclear region is interpreted as 
a jet emission (Boettcher et al. 1997). The nucleus is also 
an X-ray source with a 2-10 keV luminosity of ~ 2 x 10 40 
erg s" 1 , as confirmed repeatedly by Einstein Observatory 
(Fabbiano et al. 1988), EXOSAT (Barr et al. 1995), Ginga 
(Tsuru 1992), ROSAT (Boiler et al. 1992), BBXRT (Petre 
et al. 1993), ASCA (Ishisaki et al. 1996), and BeppoSAX 
(Pellegrini et al. 2000). These multi-mission results indicate 
a mild X-ray intensity variation. Ishisaki et al. (1996), using 
multiple observations with ASCA spanning ~2.5 years, de- 
tected long-term variation as well as a short-term (~ 1 day) 
variation on one occasion. Pellegrini et al. (2000) detected a 
periodic oscillation with a sub-day period. In addition, the 
spectra of this LLAGN, except that measured with Einstein, 
are well fitted with a power-law model having a photon index 
of 1.5-2.0. Therefore, M81 has been considered to possess an 
LLAGN. 



2 OBSERVATIONS AND DATA REDUCTION 

Owing to the explosion of the supernova SN 1993J 
(Kohmura et al. 1994) in M81, this particular field has been 
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observed with ASCA 16 times over 5.5 years, as summa- 
rized in the observation log of Table 1. It is among the ob- 
jects most frequently observed with ASCA. The results on 
SN 1993J have been reported by Kohmura (1994), Kohmura 
et al. (1994), and Uno (1997). The results on the M81 nu- 
cleus of the 1st to 10th observations have been published as 
Ishisaki et al. (1996). 

The 16 observations were carried out with the GIS 
(Ohashi et al. 1996; Makishima et al. 1996) in the PH mode, 
while the SIS (Burke et al. 1991) in several different modes 
as described in Table 1. The M81 nucleus typically gave an 
intensity of 0.3 c s _1 for the GIS (0.7-10 keV, for each detec- 
tor) and 0.4 c s -1 for the SIS (0.5-10 keV, for each detector). 
According to Kohmura et al. (1994), the SIS count rate of 
SN 1993J was ~ 25% of that of the nucleus in the 1st ob- 
servation, and decreased as t~ ' 7 , where t is the time after 
the explosion. 

For the first six observations when SN 1993J was bright, 
we accumulated photons l'.5 around the nucleus to avoid 
the contamination from SN 1993J. Since SN 1993J exhibit 
hard spectrum in this epoch (Kohmura et al. 1994), we ex- 
cluded photons above 4 keV in the following analysis. For the 
other 10 observations, we accumulated photons 3' around 
the nucleus, and considered all photons in the analysis. The 
SN 1993J contribution to the data accumulation region is 
expected to be < 5%. 

For the spectral analysis, we accumulated background 
spectra using the blank-sky data provided by NASA. The 
background photons are < 2% of the signal photons. For 
the timing analysis, we made background light curve using 
the blank regions of the M81 field. 



3 SPECTRA AND LIGHT CURVES 

Figure ^ shows long-term light curves of the M81 nucleus, 
where the 2-10 keV flux in each observation is obtained 
by spectral fitting, in the following manner. The spectra of 
the M81 nucleus are dominated by a featureless hard con- 
tinuum without significant absorption or sharp line feature, 
as Ishisaki et al. (1996) reported, so that we adopted sin- 
gle power-law model. We fixed photon index F at 1.8 and 
determined the fluxes in Figure hi, since F has remained al- 
most constant around 1.8 over 5.5 years. The GIS and SIS 
give consistent fluxes within systematic errors. The results 
on the 5th through 10th observations agree with those re- 
ported by Ishisaki et al. (1996). Through the 1st to 6th ob- 
servations conducted during the first 1.5 months, the inten- 
sity remained almost constant. On the other hand, over the 
other 10 observations performed every six months, peak-to- 
peak intensity variation amounts to a factor of three. 

We also analyzed short-term variability of the source 
intensity as shown in Figure ^, by subdividing each obser- 
vation into time bins of 1/8-day width. With the 1/8-day 
bin, typical error due to photon counting statistics becomes 
6 x 10 -13 erg cm -2 s , or ~3% of the average flux. 

We summed up the GIS and SIS counts, since the 
two instruments give consistent results. Considering the in- 
tegration radius and positional response of individual ob- 
servations and the spectral information, we converted the 
count rates into 2-10 keV fluxes. Before the conversion, 
we subtracted the background and the contamination from 



SN 1993J, using the average value of each observation. We 
fitted the light curves assuming a constant count rate, to 
evaluate the intensity variation through the \ 2 test. The ob- 
tained values of \ 2 an d degree-of-freedom (d.o.f.) are shown 
in the left bottom of the figures. Thus, the source exhibits 
statistically significant short-term variation in the 5th, 6th, 
8th, 11th, 12th, 14th and 15th observations, in addition to 
the 3rd (Serlemitsos et al. 1996) and 9th (Ishisaki et al. 1996) 
observations, as already reported. 



4 TIMING ANALYSIS 

With the 1/8-day bin width, there are 112 light-curve bins 
in total, of which the average flux becomes 2 x 10 -11 erg 
cm -2 s _1 . Root-mean-square (RMS) of the variation be- 
comes 7 x 10 -12 erg cm -2 s _1 , an order of magnitude larger 
than Poisson error. As already attempted by Ishisaki et al. 
(1996), we may estimate a gross continuum shape of the M81 
PSD using the 112-point light curve, because it includes var- 
ious sampling intervals: hours, days, months and years. 

However, the light curve is subject to an extremely 
uneven and sparse sampling. To illustrate it, if M81 had 
been observed continuously during the total observed span 
(2022 days), we would obtain 2022 x 8=16176 bins with 
the 1/8-day bin width. In practice, there are only 112 bins. 
More than 99% of the total span are thus data gaps. Conse- 
quently, it is difficult to determine the M81 PSD in general 
manner, i.e. Fourier transforming the light curve into "peri- 
odogram" , and obtaining a PSD as modulus squared of the 
periodogram. 

We have introduced two key methods to estimate the 
long-term PSD from the sparse light curve (Iyomoto 1999; 
Iyomoto et al. 1999). The first is "forward method". Namely, 
assuming various PSDs, we generate fake light curves and 
compare them with the actual M81 data to estimate which 
model PSD can best reproduce the actual M81 data. We 
utilized Monte-Carlo technique to evaluate the effects of the 
data sampling gaps. The second is the use of "structure func- 
tion", hereafter SF, which is mathematically equivalent to 
PSD, but is less affected by data gaps. We convert the light 
curves into SFs, before comparing the actual data and the 
simulation, because it is difficult to compare light curves 
directly. 

In §4.1, we calculate the M81 SF. In §4.2, assuming var- 
ious PSDs, we generate fake SFs via Monte-Carlo technique. 
In §4.3, we compare the fake and actual SFs to determine 
break frequency, /b, and slope index, a, of PSD. In §4.4, we 
determine the normalization of PSD, Po, directly from the 
RMS of the M81 light curve, and derive confidence regions 
of the M81 PSD. 

4.1 Structure Function 

The definition of first-order structure function is given by 
Simonetti et al. (1985) as 

JV-l-n 

S(rk) = -^-y ^2 ( l ™-ln+mf , (1) 

m — 

where l„ is luminosity at time t n . The sum is taken over data 
pairs satisfying t n + m —t m =Tk, with tu (fc=l,2,...) being the 
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time-lag, and N{tu) denotes the number of such pairs. For 
a stationary random process with equally spaced sampling, 
SF, auto-correlation- function (ACF), and PSD are equiv- 
alent to one another. For example, SF is related to ACF, 
A(r k ), as 



3(i%) = 2a 2 [l-A(r k )} 



(2) 



where a 2 is variance of the process. 

Nevertheless, SF is much less sensitive to data gaps 
compared to ACF and periodogram, and hence to PSD (e.g. 
Hughes et al. 1992) . In astronomy, SF have been used in var- 
ious research fields including radio (e.g. Hughes et al. 1992), 
optical (e.g. Paltani et al. 1997, Kawaguchi et al. 1998), and 
X-ray (e.g. Leighly et al. 1997; Heidt et al. 1998) frequencies. 

Figure ||a shows the M81 SFs calculated by eq. (§). In 
addition to the SF derived from the 1/8-day bin light curve, 
we superposed small-time-lag portion of SFs derived from 
1/4, 1/16, 1/32 and 1/64-day bin light curves, for compar- 
ison. The SFs with shorter light-curve bin, especially those 
with 1/32 and 1/64-day binning, are strongly contaminated 
by the Poisson noise. On the other hand, it is reasonable to 
suppose that the Poisson noise becomes negligible for the 
SFs with light-curve binning longer than 1/8 day, because 
the 1/4 and 1/8-day bin SFs are similar to one another. It 
is consistent with the fact that the Poisson noise is ~1/10 
of the RMS amplitude of the M81 variation in the 1/8-day 
bin light curve. Accordingly, we utilize the 1/8-day binning 
in this paper. 

Thanks to the various sampling intervals involved in 
our data, we can calculate the SF over a wide frequency 
range. SFs with time-lags shorter than one day come from 
within individual observations. SFs with time-lags between 
two days and 50 days are owing to various combinations 
among the 1st to the 6th observations. SFs with time-lags 
longer than six months are provided by combinations of ob- 
servations which are separated by multiples of six month. 

For time-lags shorter than one day, the M81 SF is 
smooth. Relatively large values of N(rk) may account for 
the smoothness. On the contrary, in the larger time-lag re- 
gion, the calculated SF exhibits a large scatter, due to small 
N(rk)- Therefore, hereafter, we average the SF into appro- 
priate bins in time-lag, as shown in Figure ^b. This is justifi- 
able because we are interested in the PSD continuum, rather 
than its fine structure. 

The binned M81 SF seems to keep increasing towards 
the longer end of the time coverage. In other words, it is 
likely that the characteristic time scale of the M81 nucleus 
is longer than the total observed span of 2022 days. How- 
ever, to make this inferences more secure, we need a more 
quantitative evaluation. Even the binned SF is subject to 
strong sampling bias, which depends on the sampling win- 
dow. Consequently, it is difficult to directly convert the M81 
SF into a PSD. Similarly, it is difficult to directly compare 
SFs suffering from different sampling windows. Instead, we 
here take a forward-method approach to find out a set of 
PSDs that are consistent with the observed SF. That is, we 
start from an assumed PSD, produce a number of Monte- 
Carlo light curves based on the assumed PSD and employ- 
ing the actual sampling window, and convert each fake light 
curve into a fake SF. Then, we can study how the ensemble 
of SFs derived from a single PSD behaves statistically un- 
der the particular sampling window, and whether the actual 



M81 SF can be considered as a member of this ensemble or 
not. This is carried out in the following subsections. 



4.2 Monte-Carlo Simulation 

In order to qualitatively study the characteristic time scale 
of the variation and to estimate its lower limit, we hypoth- 
esized that the M81 PSD is expressed by a convex broken 
power-law as 



P(f) 



Po 

Po (//A)- 



(/ < A) 
(/ > A) 



(3) 



as a start point of the forward method. Here / denotes the 
frequency, Po is a constant, /b is a characteristic frequency 
called "break frequency", and a is a slope index. It is an 
empirical PSD of AGNs: white at low frequency and red at 
high frequency (e.g. Pounds & McHardy 1988; Edelson & 
Nandra 1998; Hayashida et al. 1998). In this paper, we use 
this formula as a template to describe our result on M81, 
and try to estimate the values of fb, a and Po- We examine 



1.3 



whether this hypothesis is appropriate or not. We 
regard A an d ot as free parameter s of the forward method 
The value of Po is adjusted in 
1/A are 10" " 



10 J 



10 1 



4.4. The trial values of 
10 4 ' 8 days, and those of a 



are 1.10, 1.15, 1.20, 2.10. We consider Fourier periods 
from l//i = 65536 days to I//2 = 1/4 day, so that there are 
2 18 (=262144) bins in the PSD. We adopted the l//i that is 
~ 32 times longer than the total observed span of the actual 
data, considering that l//i of the actual PSD is infinite, and 
that the Fourier components longer than the total observed 
span also contribute to the observed light curve. 

Assuming a particular set of /b and a, we generate fake 
SFs in the following steps. 

(i) Using Monte-Carlo simulation, we generate a set, namely 
2 18 (=262144), of random numbers, uniformly distributed 
between and 2-n, and use them for the phase of each Fourier 
component of the assumed PSD. 

(ii) By a Fourier transform of this Monte-Carlo generated 
periodogram, we obtain a light curve of 2 19 (=524288) bins, 
with 1/8-day bin width and 65536-days total span. 

(iii) From the 524288-point light curve, we select 112 data 
points by applying the M81 sampling window. 

(iv) We normalize the 112-point light curve to have an RMS 
equal to that of the M81 light curve. 

(v) We calculate SF from the 112-point light curve. 

If we use another set of random numbers and go through 
the above steps, we can generate a different light curve from 
the same PSD. In this way, we generate 1000 fake SFs for 
each set of fb and a. 

Figure Be provides examples of fake SFs that were pro- 
duced from a single PSD having 1/A — 10 40 day and 
a — 1.4; the SFs with 1/A — 10 40 day represents those 
having 1/A enough longer than the total observed span. On 
large-time-lag range, the fake SF scatters as largely as the 
M81 SF due to the sparse sampling. Even on the time-lag 
range smaller than one day, there are significant scatters be- 
tween individual fake SFs, arising from the sparse sampling, 
while each SF is smooth in such time-lag range. These are 
explained by mutual dependences among the SFs at differ- 
ent time-lags. Namely, the smoothness in the small time-lag 
range is not only because N(rk) is large, but because there 
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are mutual dependences among the SFs at different time- 
lags. 

Using ensemble properties of the 1000 fake SFs, we can 
reduce the effect of the sparse sampling. The ensemble aver- 
age becomes smooth, as examplified in Figure |^d. Moreover, 
using standard deviation of these fake SFs, we can estimate 
how much the SF scatter under the sparse sampling. In Fig- 
ure ^d, we show the standard deviation as error bars. Thus, 
we compare each ensemble-averaged fake SF with the actual 
SF considering the standard deviation, to examine whether 
the particular PSD assumed is appropriate or not. 

It is worth pointing out that the l//i of the assumed 
PSD really have an effect on the fake SFs. Figure ^e shows 
the ensemble averages of fake SFs with l//i of 2 11 , 2 12 , 
2 18 days. The fake SFs having l//i > 2 16 days well 
reproduce "drops" around one day in the actual M81 SF. 
We confirmed that similar dependence on l//i is seen in the 
fake SFs with shorter characteristic time scale, l//b, such 
as 10 days, 100 days and 1000 days. These effect of l//i 
depends on the sampling window. Therefore, it is important 
to select 1 //1 that is enough longer than the total observed 
span of the actual light curve. In this paper, we adopted 
l//i of 2 16 day. 



4.3 Comparison between Fake and Actual Data 

In order to quantitatively compare the actual and fake SFs, 
we utilized the standard x 2 evaluation technique. We defined 

X 2 as 



2 

X = 



E 



gfakc 



(4) 



where SiPi is the M81 SF at a time-lag r k , while Sft kc 
and a^ k c are the ensemble average and standard deviation 
among the 1000 fake SFs. In this paper, we use the term 
ff^. ke t° refer to the data scattering due to the sampling 
window, and not to the Poisson statistical nature of the 
data. We regard this \ 2 as a measure of relative goodness 
of the assumed PSD. We then calculated x against various 
/b and a and obtained the minimum x 2 > namely Xmin of 
25, at l//b = 10 4 ' 6 days and a — 1.40. We can regard 
these particular set of parameters as the best estimate for 
/ b and a of the M81 PSD, as long as the PSD of eq. (|) 
is assumed. We should notice that, however, we give weight 
to the lower-limit to l//b rather than the best-fit value, 
because the obtained best-fit value is much longer than the 
total observed span, and in such l//b range, the Xmin values 
at different l//b are almost similar to one another, as seen 
in Figure [B^,. 

The next step is to estimate 68.3%, 90% and 99% con- 
fidence regions. However, we cannot apply general relation 
between confidence levels and the Ax 2 /d.o.f. values, such as 
Ax 2 /d.o.f. = 2.71 for the single parameter 90% confidence 
limits. It is because the d.o.f. is ambiguous in the present 
case, due to the mutual dependences among the SFs at dif- 
ferent time-lags. This effect is particularly prominent in the 
small-time-lag range, as mentioned in § |4.2j . 

Instead, we estimate the A\ 2 values that correspond to 
68.3%, 90% and 99% confidence boundaries, using again the 
Monte-Carlo simulation. Employing the M81 best parame- 
ters obtained above, we generated another set of 5000 fake 



light curves and hence fake SFs; we would like to emphasize 
that the Ax 2 distribution, which is given below, is almost 
independent of the assumed PSD parameters. Of course, the 
random numbers of these fake light curves are different from 
those used in § 4.2. For each of the newly produced fake SF, 



we determined the optimum /b and a, in exactly the same 
way as performed for the M81 data. We also calculated \ 2 
associated with each simulation at the assumed PSD param- 
eters themselves: l//b = 10 4 ' 6 days and a — 1.40. We denote 
it Xtrue- Note that the optimum parameters do not always 
coincide with the assumed parameters, and that Xmim cal- 
culated for the optimum parameters, is equal to or smaller 
than Xtrue- For each fake SF, we then calculated 



A 2 _ 2 _ 2 

^X — Xtrue Xm: 



(5) 



which can be regarded as a measure of statistical fluctuation 
associated with eq. (Q). Figure ^ shows the Ax 2 distribu- 
tion of the 5000 simulations. As can be seen there, 4500 out 
of 5000 simulations have Ax 2 < 23. Hence, we regard Ax 2 
= 23 as definition of the 90% confidence boundary. Simi- 
larly, Ax 2 = 3 and Ax 2 = 55 provides the 68.3% and 99% 
confidence boundary, respectively. 

Now that we have determined the confidence limits as- 
sociated with the x 2 of eq. (^), we return to the M81 SF and 
the original set of 1000 fake SFs, and examine how much er- 
rors are allowed. The two-parameter confidence ranges de- 
termined in this way is shown in Figure on the plane 
of l//b and a. The 90% confidence regions of the charac- 
teristic time scale and the PSD index are determined to be 
l//b > 800 days and a = 1.4±0.2, respectively. Our detailed 
analysis has revealed longer l//b compared to the rough es- 
timation of Ishisaki et al. (1996): they considered RMS of a 
short-term light curve (the 9th observation) and a long-term 
light curve (the 1st through 10th observations) and assumed 
a = 1.40 to obtain / b = (l-3) x 10' 6 Hz (l// b = 4-12 days). 

Now we can also examine whether the hypothesized 
shape of PSD is appropriate or not, using the Xtrue distri- 
bution of the 5000 simulations. Generally, a reasonable fit 
gives a \ 2 value similar to d.o.f. However, as mentioned in 



§ 4.2, d.o.f. is not well determined in the present case. In- 
stead, we compared the Xmin value of M81 data with the 
Xtrue distribution of fake data, as a d.o.f. independent es- 
timation of the absolute goodness. The Xtrue distribution 
of the fake light curves is shown in Figure Mb; again, the 
distribution is nearly independent of the assumed PSD pa- 
rameters of the 5000 fake SFs. Because the Xmin value of 
the actual M81 data, shown as an arrow, lies within 90% 
regions of the distribution, we conclude that the assumed 
two-component PSD is appropriate within the limitation of 
the current observation. Similarly, the single power-law PSD 
also turned out to be appropriate. 

We admit that, however, it dose not mean that we can 
exclude other shapes of PSD which we have not yet exam- 
ined. Especially, the low frequency portion of the M81 SF 
has such large scatters that other types of PSD might ex- 
plain the observed data equally well, or even better. For 
example, it may have a negative or positive index in low 
frequencies, while we assumed zero index (white noise). A 
more complex PSD, for example three-component PSD seen 
in low-mass X-ray binaries (e.g. Wijnands et al. 1999) and 
Galactic black holes (e.b. Miyamoto et al. 1994), which con- 
sists of a red-noise part at low frequency, a white-noise part 
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at intermediate frequency, and another red-noise part at 
high frequency, might explain the M81 SF as well. 



4.4 Power Spectral Density 

The remaining task is to determine Po, the PSD normal- 
ization. By definition, the variation power, J^ 2 P(f)df, is 
equivalent to RMS of the light curve with bin width of I//2 
over the total observed span of Consequently, we can 

relate Po to the RMS of the light curve, when we specify /b 
and a. 

In § 4.2, step (iv) determines Po that corresponds to 



each 112-point light curve. If there were no data gap, all light 
curves that are generated from the same /b and a would 
have the same Po. In practice, however, Po depends on the 
random numbers used in the Monte-Carlo simulation. We 
accordingly made a histogram of 1000 values of Po assuming 
l//b = 10 4 ' 6 days and a — 1.40, and determined 90% confi- 
dence region of Po. In reference to the best-estimate PQ CBt , 
the 90% confidence region in P becomes (0.8-1.2) P beat . 
Taking into account the uncertainty associated with Po, as 
well as /b and a (Figure ^|a), we have calculated a range of 
allowed PSDs with 68.3%, 90% and 99% confidences. In Fig- 
ure |^b, we plot the PSD after dividing it by the square of its 
average flux; the plot is called normalized PSD (NPSD) and 
convenient to compare random variation of different objects 
(Miyamoto et al. 1994). 

In connection with the normalization of the fake light 
curves, we would like to stress one point. In Figure ^, the 
difference in the assumed /b is best reflected in the absolute 
level of SF over the small-time-lag region; shorter values of 
l//b yield larger values of SF. However this does not mean 
that the information on longer time scales is unimportant. 
In fact, the clear separation among SFs with different PSDs 
occurs because we have normalized the RMS of each fake 
light curve. For the normalization to work properly, in turn, 
we need to know the overall data variance. When l//b is 
longer than a few times the total span, the averaged SFs with 
different values of l//b are no longer discriminated from one 
another, as shown in Figure Ha. 



5 DISCUSSION 

5.1 Comparison with Other Objects 

Let us compare the estimated PSD of the M81 nucleus 
with those of other objects. We utilize the plot invented by 
Hayashida et al. (1998): namely, a (frequency) x (NPSD) 
plot, or fP(f) plot. As Hayashida et al. (1998) pointed out, 
the fP(f) plot displays relative contribution of the power at 
each frequency decade, like the uf v plot of energy spectrum. 
As illustrated in Figure |(| the power concentrates around 
/b, so that /b is an important parameter characterized each 
system. 

So far, Ptak et al. (1998) and Awaki et al. (1999) com- 
pared the high-frequency information of LLAGNs with those 
of the other objects and suggested that LLAGNs may have 
long time scale and/or large black-hole mass. However, there 
was an implicit or explicit assumption in their analysis, that 
the individual AGNs have the same PSD shape and same to- 
tal area of NPSD on the fP(f) plot. There was no concrete 



evidence of the assumption, and it should be examined by 
measuring the low frequency portion of the PSD, as we have 
performed. Moreover, the slope index of PSD differs object 
to object (1 < a < 2). Consequently, their result strongly 
depend on the frequency where the NPSD amplitude was 
measured. Therefore, it is important to consider the PSD 
over wide frequency band. 

In Figure ^, we show the M81 NPSD, derived from Fig- 
ure |E|b, together with those of Seyfert galaxies and a Galactic 
black-hole candidate taken from Hayashida et al. (1998) and 
Pounds et al. (1988). The figure reveals two facts. First, the 
power peak of M81 locates at a very low frequency for its 
low luminosity. Second, the total area of the M81 NPSD is 
similar to those of the other objects. Consequently, the ob- 
served lack of significant short-term variability in M81 is not 
because the total variation power is small, but because the 
power concentrates on low frequencies. This suggests that 
the other LLAGNs with small NPSD amplitude at high fre- 
quency also have large power on low frequencies. 

5.2 Black-Hole Mass 

Finally, we estimate the black- hole mass, Mbh, of the M81 
nucleus. Conventionally, characteristic variation time scales 
of AGNs have been translated into Mbh assuming a linear 
proportionality between the time scale and mass, and often 
utilizing Cygnus X-l as a standard, although this estimation 
is only appropriate if the physics of accretion are similar 
between Cygnus X-l and M81. This allows us to estimate 
M BH in the M81 nucleus to be > 4.0 x 10 r M . 

So far, Mbh of the M81 nucleus has been estimated 
in various ways. From measurements of the velocity of the 
broad-line region or narrow-line region, several authors esti- 
mated the mass within these regions of the M81 nucleus, 
although their results have large scatter: 1 x 10 6 Mq- 
1.8 x 10 7 M Q (Peimbert & Torres- Peimbert 1981), (3 - 8) x 
10 5 M Q (Filippenko & Sargent 1988), and 7 x 10 5 M - 
3 x 10 6 M (Ho et al. 1996). From stellar velocity disper- 
sion measurements, Keel (1989) obtained an upper limit of 
2 x 1O 7 M0, and Bower et al. (1996) obtained a lower limit 
of 3 x 10 6 Mq. Our result falls on the relatively higher side 
of these previous estimates, though still consistent consider- 
ing large errors associated with various estimates. Our esti- 
mate reinforces the view that LLAGNs have low X-ray lumi- 
nosities because they have very low accretion rates, rather 
than because they have less massive black holes than the 
more luminous AGNs. This further encourages us to inter- 
pret LLAGNs as a relic of the past luminous AGNs. 
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Figure 1. The CIS (circle) and SIS (cross) fluxes of the M81 
nucleus in the 2-10 keV band, obtained from spectral fits with a 
photon index fixed at 1.8. Upper and lower panels correspond to 
the 1st through 6th and 7th through 16th observations, respec- 
tively. Each bin represents each observation. Systematic errors of 
the spectral fit are shown as error bars, while Poisson errors are 
negligible. 



Figure 2. The M81 short-term light curves with 1/8-day bin 
width. Abscissa is time in days after 1993 March 23 (UT), and 
ordinate is 2-10 keV flux in unit of erg s _1 cm -2 . Individual 
panels correspond to the light curves from the 1st through 16th 
observations. Observation IDs are shown in the top left of each 
panel. The vertical error bars correspond to the Poisson error. The 
horizontal error bars represent the bin width of the light curve. 
The x 2 an d d.o.f. values against the constant-count hypothesis, 
of the constant fit are shown in the left bottom of each panel. 



Figure 3. (a) The raw M81 SF without binning, made from the 
1/8-day bin light curve (crosses). For comparison, small-timc-lag 
portion of SFs made from 1/4-day (triangles), 1/16-day (circles), 
1/32-day (squares) and 1/64-day (diamonds) bin light curves are 
superposed, (b) The M81 SF after binning. The horizontal error 
bars represent the bin width of the SF. (c) Examples of individual 
fake SFs, having a = 1.40 and l// b = 10 4 ' days. Each symbol 
corresponds each SF. (d) Examples of ensemble averages of 1000 
fake SFs, having a = 1.40 and l//b=10 days (crosses), 100 days 
(circles), 1000 days (squares) and 10000 days (triangles); abscissa 
of the 100, 1000 and 10000 days SFs arc shown with slight shifts 
for convenience. The standard deviations among 1000 Monte- 
Carlo simulations arc shown as vertical error bars, (e) Ensemble 
averages of the fake SFs with various l//i of 2 11 days (circles), 
2 12 days (squares), 2 13 days (triangles), 2 14 days (diamonds), 
2 15 days (pluses), 2 16 days (asterisks), 2 17 days (crosses) and 2 1S 
days (minuses). The assumed PSD parameters are a = 1.40 and 
l//b = 10 4 6 days. The standard deviations among fake SFs are 
similar to those in (d), and are not shown for simplicity. 



Figure 4. (a) The A\ 2 distribution of 5000 fake light curves. 
Here, A X 2 = X? ruc - xLn ( b ) The X?ruc distribution of 5000 
fake light curves (histogram) and Xmin °^ ac t ua l M81 data (an 
arrow). The M81 X^in ^ cs within 90% regions of the distribution, 
confirming that the employed PSD shape is appropriate. 



Figure 5. (a) Confidence contours for the M81 light curve on 
the PSD index and break frequency plain. A cross represents the 
parameters that give the best fit to the ASCA light curve. Dot- 
ted, solid and dashed lines represent 68.3%, 90% and 99% con- 
fidence regions, respectively, (b) Estimated PSD of the M81 nu- 
cleus, shown after normalized by the average flux. Thick, medium 
and thin lines are 68.3%, 90% and 99% confidence regions of /b, 
a and Pq, respectively. 



Figure 6. Comparison between NPSD of the M81 nucleus and 
those of other objects, in the form of frequency X NPSD plot. 
Thick, medium and thin lines of the M81 NPSD are 68.3%, 90% 
and 99% confidence regions, respectively. NPSDs of Cygnus X-l, 
MGC-6-30-15 and NGC 4051 and that of NGC 5506 arc taken 
from Hayashida et al. (1998) and Pounds et al. (1988), respec- 
tively. 
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Table X. Log of ASCA observations of the M81 field. 



ID Date GIS SIS 





Start" 






End" 






exposure" 


c/s c 


exposure" 


c/s c 


Mode" 




X5 on e 


1 


93/04/05 


06 


00 


04/06 


00 


53 


32.7 


0.191 






4CCD 


B/B 


C0/C2 


2 


93/04/07 


05 


55 


04/07 


22 


12 


32.1 


0.160 






4CCD 


B/B 


C0/C2 


3 


93/04/16 


22 


40 


04/18 


17 


06 


88.8 


0.203 


59.0 


0.360 


4CCD 


F/B 


C0/C2 


1 


93/04/25 


14 


15 


04/25 


22 


05 


12.5 


0.191 


10.4 


0.346 


4CCD 


F/B 


C0/C2 


5 


93/05/01 


22 


17 


05/02 


14 


11 


24.1 


0.167 


18.7 


0.297 


4CCD 


B/B 


C0/C2 


6 


93/05/18 


20 


48 


05/19 


20 


31 


40.4 


0.210 


37.0 


0.391 


2CCD 


F/B 


C0/C2 


7 


93/10/24 


14 


33 


10/25 


14 


21 


43.5 


0.214 


29.3 


0.343 


4CCD 


F/B 


C3/C1 


8 


94/04/01 


05 


08 


04/02 


03 


21 


38.1 


0.185 


33.1 


0.327 


1CCD 


F/F 


C0/C2 


9 


94/10/21 


04 


25 


10/22 


07 


41 


48.3 


0.192 


44.0 


0.351 


1CCD 


F/F 


C2/C0 


10 


95/04/01 


18 


50 


04/02 


23 


11 


21.4 


0.281 


18.1 


0.495 


1CCD 


F/F 


C0/C2 


11 


95/10/24 


18 


40 


10/25 


22 


51 


36.4 


0.247 


35.5 


0.441 


1CCD 


F/F 


C2/C0 


12 


96/04/16 


05 


34 


04/17 


09 


47 


48.4 


0.343 


45.3 


0.541 


1CCD 


F/F 


C0/C2 


13 


96/10/27 


04 


13 


10/28 


05 


25 


32.9 


0.302 


29.5 


0.489 


1CCD 


F/F 


C2/C0 


11 


97/05/08 


12 


21 


05/09 


16 


21 


50.2 


0.315 


43.7 


0.493 


1CCD 


F/F 


C0/C2 


15 


98/04/10 


06 


18 


04/11 


10 


30 


48.8 


0.438 


38.2 


0.669 


1CCD 


F/F 


C0/C2 


16 


98/10/20 


19 


54 


10/21 


23 


11 


5.8 


0.240 


15.2 


0.334 


1CCD 


F/F 


C2/C0 



a The start time of observation in year/month/day 
houriminutes, and the end time in month/day hounminites. 

b Exposure time in seconds after the standard data screening. 

c Observed count rate of the nucleus in 0.5—10 keV. 

d The SIS data mode for high/medium bit rates, with F and 
B denoting Faint and Bright mode, respectively. 

e The chips of the SIS0/SIS1 detectors on which the M81 nu- 
cleus was placed. 
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